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1. Introduction 


In recent years increasingly complex ving planforms, which are 
0 f“ficient over a large flight envelope , are being used in the aircraft 
industry. An efficient high-subsonic or supersonic cruise aircraft 
must be designed to work efficiently even for off-design performance 
points, such as, landing and take-off. One way of designing such an 
aircraft is to have fully attached flow at cruise and controlled 
leading— edge separation at landing— and take— off— conditions (ref. 1 ). 

A significant amount of vortex lift can be generated at high angles 
of attack by leading-edge vortex flow. Henderson (ref. l) pointed 
out that at high angles of attack highly swept-back wings having sharp 
leading-edges have low amount of leading-edge suction (essentially 
zero) and generate large amount of vortex lift. However, the 
leading-edge suction increases as the sweep is decreased or leading- 
edge radii are increased. The structural considerations restrict the 
leading-edges to be of finite radii. In these cases the lift data 
lies between potential and potential plus vortex estimates (zero 
suction). At present, there exists no theoretical method which can 
predict aerodynamic characteristics of wings having partial leading- 
edge suction and vortex flow as mentioned above. A theoretical method 
is presented here to predict aerodynamic characteristics of wings 
under such conditions. 
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A theoretical method for cases vith non-zero leading-edge suction 
and partial leading-edge separation should also he applicable to cases 
with complete leading-edge separation. Due to lack of data vith partial 
leading-edge separation, the present method has been extensively 
compared vith other theoretical methods and available data only for 
cases with complete leading-edge separation. Thus, a survey of 
literature on completely separated leading-edge vortex flow seems 
pertinent . 

Legendre (ref. 2) va? the first one who attempted to solve the 
leading-edge separation problem. He assumed that rolled up vortex 
sheets over the wing can be replaced by a pair of concentrated vortices 
(fig. l). In this model the flow tangency boundary condition over the 
wing w'as satisfied along vith the leading-edge Kutta condition. The 
flow was assumed to be conical so that conformal mapping could be used. 
Upon the instigation by Adams (ref. 3), Legendre (ref. revised his 
model. He reported that in his first model he had assumed a cut 
Joining the two vortices. In another form of his model, Legendre 
included a cut between the vortices and their respective leading-edges 
tc account for the feeding sheets. Adams pointed out that in the 
first nicdel cf Legendre the lift was multivalued because the region 
was no longer simply ccnnected and in the second one the pressure 
difference was allowed across the sheet. Based on the suggestions of 
Edvards (ref. 5 ), Brown and Michael (ref. 6) modified Legendre's 
slender tody model by using feeding cuts, which ccnnected the line 
vortices to the wing leading-edges (fig. 2). This vortex system of 



3 


concentrated line vortex and feeding cuts was required to satisfy 
the force free condition. The leading-edge Kutta condition and 
flow tangency boundary condition on the ving were also satisfied. 
Mangier and Smith (ref. 7) proposed a somewhat more realistic model 
than that of Ercvn and Michael, but still used slender tody theory 
in their investigation. They used a continuous model of the separated 
vertex sheet along with a concentrated core (fig. 3). The shape emd 
strength of the vortex sheet and the concentrated line vortex were 
determined by satisfying the flow tangency condition on the ving and 
the pressure continuity condition across the separated vortex sheet. 
Later, Smith (ref. 8) used segmented feeding vertex sheet with 
considerable improvement in numerical procedure (fig. ^), largely due 
zc advent of greater computing power. The above model vas modified 
further for thick vings (Smith, ref. 9). The main shortcoming of all 
these models described sc far is the assumption of conical flow. 

Gersten (ref. 10) extended Bolla^^’s vortex model (ref. 11), which 
was for rectangular wings with ving tip separation, to arbitrarily 
shaped vings of small aspect ratio with leading-edge separation 
(fig. 5). In this model the vortices came off wing edges at an angle 
ci/2 zo the ving plane. The ving was replaced by infinitesimal lifting 
elements and the strength of vortices was assumed to var^^ along the 
span. At this point the flow tangency condition on the ving was 
satisfied to find the wing characteristics* Garner and Lehrian (ref. 
12) followed Gersten ’s approach by using Multbopp’s lifting 
surface theory (ref. 13) to represent the ving. Both of these models 



k 


are very crude and give only total characteristics of the ving; i.e., 
pressure distributions are not ceilculated. 

Sacks, et al (ref. 1^; ) assumed aerodynamic characteristics (lift 
and pitching moment) to be composed of two components - the lineeir 
and non-linear components. The linear component was calculated by 
using the integral method of Lawrence (ref. 15), while the non-linear 
component was calculated by assuming that vortex pairs were shed just 
outside the wing leading-edge with shedding rate determined by either 
an empirical method or slender wing theory. The location of each 
vertex pair was determined by satisfying the force free condition at 
the vortices. 

Nangia and Hancock (ref. 16) extended Brown emd Michael's model 
(ref. 6) to non-slender wings. In their model the wing planform was 
represented by bound- and trailing-vorticity distributions. The 
wake behind the trailing-edge was free to move outboard and the 
leading-edge separation was represented by two isolated vortices, 
which were connected to the leading-edge by cuts (as in Brown and 
Michael). The Kutta condition was satisfied along the leading- and 
trailing-edges and the flow tangency condition was satisfied on the 
wing surface. Zero force condition was satisfied on the isolated 
vortices and the cuts at selected collocation points. Although this* 
method was not restricted to slender wings, the leading-edge flow 


model- was crude. 
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Polhamus (ref. IT and l8) used leading-edge suction analogy to 
predict lift coefficients for various simple planforms, such as > 
arrow- , diamond- and delta-wings. It was essiimed that when the 
complete flow reattachment occurs inboard of the leading-edge vortices, 
the total lift equals the sum of the potential and vortex lift. 

These components of lift were calculated by using a modified form of 
Kulthopp's lifting surface theory (ref. 19)* The vortex lift was 
assumed to be equal in magnitude to the potential flow leading-edge 
suction force lost due to separation. In its original, form this method 
did not calculate the local distribution of lift and so the pitching 
moments were not predicted. Snyder and Lameu: (ref. 20) used this 
method to predict the longitudinal load distribution and pitching 
moment for delta wings. 

Mook and Maddox (ref. 21) modeled the leading-edge vortex system 
by finite vortex elements coming off the leading-edge (fig. 6). This 
network of vortex elements was superimposed on the vortex -lattice 
used by Giesing, et al (ref. 22). The solution is obtained in an 
iterative manner by satisfying the flow tangency boundary condition on 
the wing surface, approximately satisfying leading- and trailing-edge 
Kutta conditions, and satisfying force free conditions on the vortex 
elements over the wing surface. The force free condition was not 
satisfied on the veike behind the trailing-edge. Kandil, et al (ref, 23) 
modeled the flow in a manner similar to that of Mook and Maddox 


(ref. 21) and extended it to vlng-tip separation also. Kandil, et al 
followed Belotserkovskii (ref. 2k) for the representation of the 
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wing surface. In this model the bound elements of the vortex-lattice 
were unsvept , and the wake behind the trailing-edge and the vortex 
elements coining from wing-tip were force free. Apparently, this 
method has beer, restricted tc an angle of attack of 20 degrees (ref. 

25 ). Kandil, et al (ref. 26) extended their model to calculate the 
location and strength of a concentrated cere, which they also used 
for convergence criteria. Rehbach (ref. 27) also followed 
Eelotserkovskii to model wing-tip separation and the approach was 
similar to that of Kandil, et al (ref. 25). However, he solved the 
leading-edge separation problem differently (ref. 28). The process 
was started by finding a converged solution for a rectangular wing 

(fig. 7). The leading-edge span of the rectangular wing was decreased bv 
a small amount, while the trailing-edge span was kept constant. A new 
converged solution for this wing vas ettained. This process vas 

repeated until the planform reduced to e delta ving. The deficiency 
of these methods is that the leading-edge Kutta ccndiiion is only 


epprcximately satisfied. The iteration process of Rehhach could also 
be quite time consuming. 

Nathnen (ref. 29) presented tve models of leading-edge separation; 
tic fixed nodel and the free vake r.odel (fig. S). In oozh models 

the wing was represented by panels with constant strength doublet distri- 
bution, being equivalent to closed jvortex filaments on the boundary of 
the panels. In the fixed vake model the separated sheet war modeleti bv 
placing a ?f planar boxes along the leading^edge which extended 
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to the vortex core predicted analytically hy Brovr and Michael (ref. 6j. 
The doublet strengths were obtained by satisfying flow tangency 
boundary conditions only. In the free v&Xe model the tieparated 
vertex sheet vas represented by discrete vortices attached to the 
leading-edge. These vortices were aligned along the local velocity 
vector to be force free. The valce behind the trailing— edge was force 
free only to a certain extent. The fixed wake model is too crude, 
whereas no definite convergence criteria has been established for 
the free wake model. 

The most sophisticated and realistic model of all leading-edge 
separation models has been due to Brune, et al (rel. 3C). In this 
model the wing and the separated vortex sheet were represented by 
piecewise continuous doublet distributions. The separated vertex 
sheet vas connected to a concentrated cere by a cent mucus fed sheet 
(fig. y)* The solution vas obtained in an iterative manner b^- 
satisficing the Kutta condition along all edges, flow tangency boundary 
conditiori on the wing, and the force free condition on the sepa.rated 
vertex sheet and vake behind the trailing-edge. This model has also 
teen extended tc thick and cambered wings. The drawbacks of this 
medei are that it cannot predict lift correctly at small OTiglcs of 
attack for moderate to low aspect ratio vings , takes toe much 
computer time to get a converged solution and needs large computer 
ceDory space . 
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All methods mentioned above celcLaate only the wing characteris- 
tics for complete loss of leading-edge suction. A method is 
developed here for partial leading-edge flow separation with non-zero 
leading-edge suction. The wing boundary condition is formulated by 
the Quasi Vortex Lattice method (QVLM) of Lan (ref. 31). The advantage 
of this method is that the leading-edge boundary condition can be 
exactly satisfied. Tlie leading-edge separated vortices are 
represented by discrete free vortex elements which ore aligned with the 
local velocity vector at their mid-points to satisfy the force free 
condition. The vahe behind the trailing-edge is also force free. The 
flow tangency boundary condition is satisfied on the wing, including 
the leading- and trai ling-edges. Due to the non-linear nature of the 
problem, the problem is solved in an iterative manner. Eue to non- 
availability of any data with partial leading-edge separation, the 
method will be compared only with other methods (ref. 32 and 33) and 
experimental data (ref. thru ko) for complete leading-edge 
separation. The basic assumption in the present theory is that the 
Prandtl-Glauert equation is applicable. The thickness and fuselage 
effects are ignored. 

Chapter 2 presents the theoretical method. In Chapter , 
numerical results are presented and discussed. Conclusions and 
recommendations are made in Chapter U. 



2. Theoretical Method 


2.1 Problem Definition 

In steady symmetric flight at a high angle of attack, the flov 
over a thin low aspect ratio highly sveptback ving separates along 
the leading-edge and the tips. In the folloving, only delta vings 
will be considered. The ving can be represented by a bound vortex 
sheet, across which exists a pressure difference, and the separated 
flov along leading-edges by force free vortex sheets, across which 
there is no pressure difference. In the present method, the Quasi- 
Vortex-Lattice method (ref. 31 ) is used to simplify the induced 
velocity expressions due to the bound vortex sheet and discrete force 
free vortex elements for separated vortex sheets. 

The folloving boundary conditions are imposed on the flov model; 

a. The flov must be tangential to the ving camber surface. 

b. The leading-edge boundary condition and trailing-edge Kutta 
condition are to be satisfied. 

c. The vortex elements over the ving and vake behind the 
trailing-edge are force free. 

This is a non-linear problem because the strengths cf the ving hound 
vortices oiid free vortices, and the locations of the free vortex 
elements are urltncwn. Thus, the problem is solved by an iterative 
method. 
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2,2 Mcdel Geometry 

The origin of the rectangular coordinate systeir. is at the ving 
apex. ‘The ving lies in the x-y plane and the x-axis is tahen along the 
ving center-line. The ving span is given by h and the surface area 
by £. 

2.2.1 V/ing Geometry 

The location of bound- emd trailing-vortex elements for e typical 
case are shovn in figure ]0, a detailed description of which is given 
in section 2.3. The x-location of bound elements is given by the 
cosine lav and is illustrated in figure 10. 


X . 
1 


+ 2. r- 


- Ccs( P - - ^ r)) 


( 2 . 1 ) 


1 , 2 , — 


where is the leading-edge x-coordinate, c is the chord and N is 

the number of bcurid elements in a chcrcvise direction. The spanvise 
location of trailing elements is given by, 




tt)) 


( 2 . 2 ) 


J = 1,2, H 
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where h is the span and M is the number of legs of trailing 
vorticity, which is one higher than the number of spanwise strips of 
bound elements. The locations of control points are given by, 

j 

k = 0,1,2, N 

^cp = r 
J 

J = 1,2, (M - 1) 


where 



X p and c . 
)L. 0 

respectively. 


are the leading-edge x-coordinate and chord at 


2.2.2 Leading-Edge Vortex System Geometry 

The leading-edge vortex system is superimposed on the regular 
quasi -vortex-lattice grid. A typical vortex element is shown by 
points A through J in figure 11. These. points are connected by a 
series of short straight segments. The initial location of these 
segments is shown by dashed lines and final location by solid lines. 
These segments have the following characteristics: 


OF POOR QUALOT 
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a. Points A through E lie along a ving trailing vortex element . 
Initially point A is one root chord, away from the trailing- 
edge in the downstream direction and the line segments 
"between A and D are parallel to the axis of sjTnmetry. The 
line segments "between points A and B are of equal length. 

In the final converged position these segments are aligned 
in the direction cf the local velocity vector. The segments 
B-C and C-D are 0.1 Cp long. B-C is allowed to move only 
in the vertical direction whereas C-D is fixed in the wing 
pleuie because flow is tangential tc the trailing-edge. 

Segment I>-L is also fixed in the ving plane. 

b. Points E, F, G and K also lie in the wing plane. The location 
of segment E-F is ahead of the ving first bound element and 

is given by. 


X 


E 


X 


F 



+ ~ (1 - Cos( 
~ (1 - Cos( 

C 


- - " \] 

2(N + 1)^^ 

I—)) 

2(N + 1)^' 


(2.5a) 

(2.5b) 


where the subscripts E and F refer to the points under 
consideration. The above two equations are similar to 
equation (2.1), It is to be noted that segment E-F is 
located at the first bound element for a grid of (N 1) 
bound elements in a chordvise direction. The segments F-G 
and G-H are of the same length and point G lies on the 
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leading-edge. The segment G-H is fixed in the wing plane 
due to the leading-edge boundary condition, 
c. The initial location-of point I is given by, 


X 


I 


(2.6a) 


“ y . ( 2 . 6t ) 

I '’I'’ 

7 .^ = 0.1 C ,3 tan(22.5 - 0.5a) for a < 15° (2.6c) 

I n 

or z., = 0.1 C_ tan a for a > 15° (2.6d) 

I r. 


where C„ is the root chord and a is the angle of attack. 

n 

Initially point J is one root chord away from the trailing-edge 
The segments between point I end J are of equal length and 
lie in a plane parallel to x— z plane. These sepnents are 
approximately at a height of 0.1 Cp above the wing plane (see 
Chapter 2.6). In the final converged position all the segments 
between points H and J are aligned in the direction of the 
local velocity vector. 

d. The semi-infinite segments from points A to infinity and J tc 
infinity are straight and are parallel to the undisturbed 
free-stream direction. 
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2.3 Induced Velocity Due tc Wing 

In the quasi-vortex-lattice method (ref. 31 ), for the purpose of 
satisfying ving "boundary (tangency) condition, the continuous vortex 
distribution over the ving is replaced by a quasi-continuous one, 
being continuous chordwise but stepwise constant in the spanwise 
direction. Thus, the wing surface can be divided into, a number cf 
vortex strips with the associated trailing vortices (fig. 1 C). Ir. an: 
strip, consider a vortex element ydx with an arbitrary direction £ 
(fig. 12 ). The induced velocity due tc all bound elements in ith 
strip is given by (see Appendix A), 


,2 x 


I- (R) = Y(x’)f— 


t 


X 


-7^} • £^dx' ( 2 .T) 

la- I I 


and due tc the associated trailing vertices by (ref. 31), 


.2 X, 


-*■ R‘- - R) X dl 

= fc C. •'x. V 1 


( 2 . 8 } 


c ( V ) 

The transformation, x' = x, + - y -a (i - Cos 6), reduces ecuations 


( 2 • 7 ) and (2.8) to, 


q_ (R) = G^( 0 ) Y(e) Sin 6 d6 (2.9) 

1 


and 


\ (R) 

2 


^ Aiy) 

Bv 'o 


02(6) y( 8) Sin 0 d 0 


(2.10) 
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where , 


G^(e) = 


a X 


la' X i.' 


1 _ iJl 


G^(6) = 


(R. - R) X d£ 


t' 


and c(y) = \ ~ affined in the List uf 

The total induced velocity due to the ith strip of vortex distribution 

is given by , 


q^(R) = G^(e) y(6) Sin 9 d9 

2 

+ G!(e) y( 6) Sin 6 d9 

8tt o 2 

_ ^7T g.^ e (2.11) 

OTT O 2 


where the first tern is due to bound elements, second due to left leg 
of trailing elements and third due tc the right leg of trailing elements. 
The above integrals are reduced to finite sums through the midpoint 
trapezoidal rule (see ref. 31 )t 


u (B) > i: <5, ♦ 2; - 2; sin e. 




(2.12) 


where 


0 _ ( 2 k - .Ij. ^ locations of boxmd elements are given by. 


2N 
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k 1 


= |- [1 - tt)], k = 1,2, ,N 

X. = the leading edge x-coordinate at (left leg) 

^1 ^ 

X. = the leading edge x-ccordinate at y (right leg) 

2 ^ 

= chord length at y^^ 

Cg = chord length at y^ 

The control points in the chordvise direction are chosen such that. 


(2.13a) 

(2.13t) 

(2.13c) 


X . 

1 



-os — ). 1 


1,2, N 


(2.14) 


where x^ is the leading edge x-coordinate on the chord c through the 
control point. The spanvise location of trailing vortices is given by, 


y . 

j 



Cos i 


- 1 


j , 




— M 


(2.15) 


and control points by, 




= 1,2, (M-1) 


(2.l6) 


where b is the span Eind M is the total nuinber of trailing vortices 
which is one more than spanvise strips. The geometry associated vlth 
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equations (2.13) - (2.l6) is based on the semi-circle method and is 
illustrated in figure 10. 

Thus, induced velocity due to all vortex strips of the wing can be 
written as , 

M-1 

q(R) = E q.(R) (2.17) 

i=l ^ 


2 . h Induced Velocity Due to Leading-Edge Vortex System 

The leading-edge vortex systeri, as described in Chapter 2*2.2, 
consists of elements. Each element ma\' have different number of 

small vortex segments. Assume that ith set has L small segments. 

The induced velocity at a point (x,y,z) due to Jth segment of ith 
element is given by (Appendix A), 


V. (E) = 


2 

6 r. 

2 

hv 




► "t 1 2 I 

l' X V I |b' I la’ 


. t’ 


( 2 . 18 ) 


where 


R = xi + + zk 


a = (x. - x)i + (y. - y)J + (z, - z)k 
J J 


I = - x)l > (y,^^ - y!j ^ - Xik 
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a' = (Xj - x)i + B(yj - y)j + 6(z^ - z)k 

b’ = (x^^^ - x)i + - y)J + - z)k 

!■ . - y.)J * 

= Vortex strength of ith set of segments 

The subscripts J and j+1 correspond to the end points of Jth 
segment . 

Now the induced velocity due to ith element can be written as, 

V (R) = I V (R) (2.19) 

^=1 j 

Therefore, the induced velocity due to all elements is, 

^ ^ M-1 ^ 

V(R) = Z V (R) (2.20) 

i=l 

2. 5 Boxmdary Conditions 

The two basic boundary conditions to be satisfied in the model are, 

a. The flow must be tangentiaJ to the wing camber surface. 

b. The vortex elements above the wing and in the wake behind 
the trailing-edge must be force free. 
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2. 5*1 Formulation Of Ving Boimdary Conditiop 

The bound elements and the corresponding control points of the 
ving surface are numbered from the leading-edge to the trailing-edge 
and ft-om the root to the tip. Thus, there are Na = N(M - 1) bound 
elements and corresponding control points (see Chapter 2.2.1). 
Similarly, tne vortex elements and the corresponding control points 
in the leading-edge vortex system are numbered from the root to the 
tip. There are Ifb = (M - l) leading-edge vortex elements being 
equal to the number of wing vortex strips. Thus, there must be 
Nc = Na irb control points on the wing siirface. The flow tangency 
condition can be written as. 


r 








NcxNa 

Waxl 

NcxNb 

Nhxl 


(||) - Sin a 

= J 




k k 

Nc X 1 


J 



( 2 . 21 ) 


where A , is the induced downwash at ith control point of wing 
ij 

due to a unit horseshoe vortex density at the induced 

downwash at ith control point due to kth leading-edge vortex element 

of unit strength; y the vortex density of the Jth hound element; 

J . 

the strength of kth leading-edge vortex element, the camber slope 


at the ith control jwint and a the angle of attack. According to 


equations (52) and (53) of ref. 31, the leading-edge thrust coefficient 
is related to the nonnalvash on the leading edge by the following 
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equation : 


a. = induced ncrmalvash - (§^) * Sin a 

k 


( 2 . 22 ) 


which leads to the last expression on the right hand side of equation 
( 2 . 21 ). In the above expression, is defined as: 



( 2 . 23 ) 


N = Number of chordvise vertex elements 
A = Leading-edge sweep angle 
= Free stream Mach number 

CO 



= Sectional leading-edge thrust coefficient of the kth strip 

For complete leading-edge separation cases the sectional leading-edge 

suction is zero and sc is . 

"k 

2.5*2 Fonnulation Of Force Free Condition Of Free Elements 

The vortex segments above the wing s'orface and the vajLe are to be 

aligned in the direction of local velocity vector calculated at their 

mid-points. Consider ith segment of a vortex element. The 

coordinates cf its end points are given by (x^. a y. , z.) and 

til 

^i+1’ Assume that the velocity at the mid-pcint of this 


segment at a given iterative step is given by. 



Then, the nev location of the (i + l)th end point vill be. 


(2.2Ua) 

(2.2Ub) 

(2.2Uc) 


where 

and 


-p 


2 P 2 

U = ,/u 4 + V 


As 


= i/(Vi " ^i^^ ^^i+1 - ^i'^ ( 


^-.1 - 


^-^1. 


U A 

Xi + - As 


^i+1 


N 


yj j As 


"i-i 


N 


Z, 4 — As 
1 L 


Before equations (2.21) are used, the following points sho\ild be 
considered: 

a. The. length of each segment is to be preserved. 

b. freo vortex segments above the wing sho^jJLd net come toe 

close to the wing surface to avoid numerical difficulty 
in the present inviscid theory. 

c. The adjustment of the location of each segment to satisfy the 
force free condition should be such that it does not result 
in numerical fluctuations (see Chapter 2.6). 




IS 
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Based on the above consideration, equations (2.2b) viil be modified 
as follows. 

Consider the same ith segment. If this segment meves a-percent 
only according to the velocity computed at its mid-point, then 
equations (2.2^) can be modified to be. 


/ly,. = ~ Ls + (1 - a)(y. , - y.) (2.25a) 

o 1 

av 

" “u " ^i^ (2.25b) 

- Ay^~- ( 2 . 25 c) 

It follows that , 


' ''i * 


(£.£6a) 



^i " 




(2.26b) 



z. + Az_ 
1 W 


( 2 . 26c ) 


Let minimum vertical distance ary vortex segment is 

allowed to come close to the wing surface. If 2 ^ . is less than 

^t is then set equal to Z . and Az is recalculated by using, 
mm min IN ^ » 


Az 


N 


Z , - z. 

nin 1 


(2.27) 
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This veilue of Az„ is used to calculate x. ,, • 

N 

2.6 Solution Procedure 

The hasic unknowns of the prohlem are the bound vortex density 
on the wing, and the strengths and the locations of the elements of 
the leading-edge vortex system and the wake. The problem is non- 
linear because the locations of the leading-edge vortex system and 
the wake are unknown a priori. Therefore, the problem will be solved 
by the iterative process described below; 

a. Prescribe the vortex lattice for the wing surface, and the 
initial locations of the free elements over the wing and 
in the wake. 

b. By satisfying the wing boundary condition, i.e. equation 
(2.21), obtain the bound vortex density of the wing and the 
strengths of free elements. 

c. Calculate all the aerodynamic characteristics, 

d. Calculate the forces acting on the free elements over the 
wing surface. 

e. Adjust the free elements of the leading-edge vortex system 
and the wake in the local velocity vector direction, as 
described in Chapter 2.5.2. 

f. Repeat steps b throiigh e until a converged solution is 
obtained. 

The initial locations of the free vortex elements are assumed by 
letting them leave the leading-edge in the undisturbed free-stream 
direction up-to a height of about ter. percent of the root chord 

PAGE 

QUil T TTV 


ORIGINAL 
OF POOR 



2k 


beyond which the elements are parallel to the wing plane. Initially, 
all the elements of the wake lie in the plane of the wing. In the 
iteration process, the force free condition is satisfied on the free 
elements from the root to the tip in the down-stream direction. A 
similar approach has been used by Butler and Hancock (ref. i<l) with 
success for the wake problem. The elements over the wing are adjusted 
before the elements of the wake. In the first iteration the segments 
over the wing are moved ICC percent according to the velocity computed 
at their mid-points. This movement is gradually reduced in steps of 
90, 80 and 75 percent in the next three iterations, after which it 
remains at 75 percent (see equations (2.25) and (2.26)). The segments 
in the wake are moved only 50 percent in each iteration. Thus, exact 
force free condition is not enforced because whenever ihe free elements 
cone close to each other they induce unreasonably large velocities 
because viscous effects are not included in the present theory. These 
large velocities increase the forces cn the elements and induce 
fluctuations in their locations. 

The solution is assumed to have converged if in two consecutive 
iterations the difference between the total strengths of leading-edge 
free vortex elements is less than one percent ana the absolute force 
acting on the free elements is in the neighborhood of e minimum. Thus, 
an exact force free condition is not enforced as discussed in the 
previous paragraph. 
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2.7 Aerodynamic Characteristics 

The expressions for the evaluation of the pressure distribution 
are derived in Appendix B. They are obtained by applying the Kutta- 
Jovikowski theorem to the vortex system on the wing. 

The sectional normal force coefficient of Jth strip is given by 
chordwise integration of the differential pressure coefficient: 

^t 

c = — / ■' iC 4x (a. 28) 

"j "j ‘n ^ 

J 

where x_ and x are the leading and trailing edge x-coordinates 

J 

of the chord passing through the control points of the Jth strip and 

c. is the chord length. The transformation, 

J 


X 




- Cos 6) 


(2.29) 


reduces equation (2.28) to; 


C = 4 AC Sin e d6 
. o F 


N+1 


-7- rr ^ AC Sin e. 

k=l *■ 


( 2 . 30 ) 


and 
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where the midpoint trapezoidal rule has been used to reduce the 
integral to a finite sim. Similarly, the sectional pitching moment 
coefficient for the jth strip about the y-axis is given by. 


' 1 

; r / ^ iC X dx 


N+1 


^ ^ I AC (Xp + ^ (1 - Cos e })Sin 6, (2.31) 

2C(N + 1) k=l ^ ^ 


where C is the mean geometric chord. 

According to equations (52) and (53) of ref. 31, the sectional 
leading edge thrust coefficient is given ty: 


Jl - Cos^A (w! - - Sin a))' 

Jl ^ .1 dx .1 

2N^ Cos A(1 - + tan^A) 


( 2 . 32 ) 


where is the free stream Kach number, A tlie sweep angle of the 

d 2 

leading-edge, and v’ and (~) are the induced normeLlvash and slcne 

j j 

of the wing surface at the leading-edge. 

The normal force coefficient is obtained by integrating the 
sectional ncrmal force coefficient across the span: 


b 

Sj ' K' '=n‘= 

J ^ 

“ 2 


(2.33) 


where b is the span and S the wing area. By the transformation, 
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y = ^ (1 - Cos <p) 


(2.3i+) 


equation (B.21) can be reduced to. 


S ‘ ^ ^0 'n ' 


M-1 

= I c c. Sin 0. 

2SM . , n. 1 2 

1=1 1 


4-i = I IT. i = 1,2, (M-1) 


U. ) 


where (M-1 ) is the total number of spanwise strips and the regular 
trapezoidal riile has been used. 

Similarly the pitching moment and leaxiing-edge thjnist coefficients 
are given by; 
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^ M-1 

bTT _ 

2SM . , 
1=1 


bTT 


M-1 


2SM . - 
1=1 




Sin 

Sin 


(2.36) 


^ , ) 


The normal force coefficient and leading-edge thrust coefficient 
can be resolved in the free stream direction and perpendicular to it 
as shewn in figure 13 to obtain the lift coefficient and the induced 


drag coefficient: 


CD 


C, = C.. Cos G 
L N 

a 

(2.38a) 

Sin a 

u w 

a 

(2.38b) 

Sin a 
a 

(2.38c) 

“ C|-j-i Ccs G 
t 

V iL , JC-U , 

where a is the angle of attack, Equetions (2.38) 

cem now be used to 

obtain the total lift and induced drag coefficients: 


C, = C Cos G + Cr- Sin a 

Lj h 

(2.39) 

Sin a - Cos a 

i 

( 2 . it C y 

in'hen the flow along the leading-edge is completely s 

eparated, the 

leadiiig-edge thrust coefficient is zero. 




3. I^esults and Discussion 


It has been found during the investigation that the calculated 
induced velocities due to the ving become inaccurate if the control 
point of a free vortex segment, \diere induced velocities are to be 
evaluated, is any closer to the ving than twenty percent of the local 
chord. On the other hand, the induced velocities calculated above the 
control points of the ving show a smooth trend. Therefore, if a free 
vortex segment is closer than twenty percent of the local chord, the 
induced velocities at its control point, i.e. mid point due to the 
ving, are obtained' by linear interpolation of the velocities calculated 
above four ving control points among which the point is located. 

It has also been found numerically that the aerod^sT.er.ic 
characteristics depended on the number of spanvise strips, i.e. M of 
equation (2.2). Therefore, a parametric study has been made to find 
a relation between the aspect ratio and the number of spanvise strips 
for reasonably accurate results (Fig. 1^). It is to be noted that 
as the aspect ratio is decreased, the number of spanvise strips has 
to be increased. This is due to the fact that the spanvise variation 
of aerodynamic characteristics, such as pressure coefficient and thrust 
coefficient, is large for small aspect ratio vings. This study was 
performed by matching the lift coefficients obtained by using the 
present method to those obtained by using suction analogy (ref. 32) 
at one angle of attack. 
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The free elements of the leading-edge vortex system heve been 
restricted not to come any closer than a minimum specified height to 
the wing surface, which is given empirically by, 

Z = 0.1 C tan(22.5 - 0.5a) for a < 15° (2.bla) 

for a ^ 15° (2.iilb) 

where Cp is the root chord and a the angle of attack. This 

restriction was needed because whenever the free elements are close to 

the wing surface, they induce large velocities on the wing and vice 

versa, which makes the free elements fluctuate (unstable). In the 

real flow, at small angles of attack, the leading— edge vortex system 

is weak and diffused. The present method does not account for 

diffused vortices and so the effect of the free vortices is artificaliy 

reduced by increasing Z . as the angle of attack is decreased below 

min 

15 degrees. 

All the results have teen calculated by using six chordvise 

vortex elCTents on the vine; i.e., N of equation (2.1) is 6, and the 

length of the free vortex segments being 15 percent of the root chord. 

The effect of the number of chordvise vortex elements and the length 

of the free vortex segments is insignificant. 

A computer program has been developed for the present model vith 

the above restrictions (ref. ^2). It has been used tc generate 

aerodynamic characteristics for flat delta vings of several aspect 
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Z . = O.I tan a 

mm E 
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ratios. These predicted results are compared vlth the avaij.atle 
experimental data and the results by the suction analog^" (ref. 32), 
Kandil’s model (ref. 23) and Brune’s model (ref. 30) as obtained by 
Kuhlman (ref. 33). The lift- and pitching moment-coefficients are 
plotted against angle of attack for complete leading-edge separaticn 
cases, i.e, zero leading-edge thrust, in figures 15 through 19- In 
general the agreement for the lift coefficient betveen the present 
method, suction analogy (ref. 32), Brune's model (ref. 30) and 
experimental data is quite good. The present method usually over- 
predicts the lift coefficient at small angles of attack vhereas 
Bruners model (ref. 30) underpredicts it. For the wing of aspect 
ratio 0.7053, the present irethcd becomes less accurate at high angles 
of attack (figure 15). This could be due to the large rate of change of 
pressure coefficients in spanvise directions at large angles of attack 
for small aspect ratio vings. A better agreement could be obtained 
by increasing the number of spanvise strips for small aspect ratio 
wings at large angles of attack. Pji excellent agreement is seen fcr 
the pitching moment coefficients calculated by using the present 
method and the experimental data in figures l6 and IT* The suction 
analogy can not predict accurate pitching moment coefficient because 
it does not calculate the surface load distribution. Although 
surface load distribution is predicted in Brune's model (ref. 30), 
the pitching moment coefficients are not predicted accurately . The 
pitching moment coefficients predicted by Kandil's model (ref. 25) for 
aspect ratio 2.0 wing are in a better agrcciucnt with experimental data 
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than any other method, hut the model is restricted up to a 20 degree 
angle of attack only. The effect of Mach number on the lift- and 
pitching moment coefficients at different angles of attack for aspect 
ratio 1.5 ving is shown in figure 20 eind the trend for the lift 
coefficient agrees with that predicted by the suction analog^'. 

The pressure distribution for three delta wings at several angles 
of attack and constant x— locations are shown in figures 21 through 23. 
In general the pressure peak obtained by using the present method is 
lower than the experimental value and is shifted towards the root 
chord. Figure 22 shows the only comparison with the theoretical method 
of Brune (ref. 33) and a sharp peak is visible in Brune's moael. 

The reason for the peak being lower in the present model is that each 
free vortex acts as a concentrated core by itself whereas Brune 's 
model has a separated vortex sheet w^ith a concentrated core at its end. 
Therefore, in Brune 's model a sharp pressure peak will be present, 
whereas in the present model the pressure distribution will be more 
diffused. 

Thus far it has been shown that the present model gives reasonable 
results for completely separated flew along the leading-edge. The 
theoretical effect of partial leading-edge separation on the aero- 
dynamic characteristics will be shown next. Figure 2k shows the 
effect of varying the amount of the leading-edge suction lost on the 
aerodynamic characteristics for delta wing of aspect ratio 2. It can 
be seen from the figures the.L fer a fixed single of attack, the lift 
coefficient and the induced drag increases as the amount of the 
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leading-edge suction lost increases. These trends are similar to the 
ones shown by Henderson (ref. 1). For these cases the lift is found to 
be highly nonlinear with angle of attack. On the other band, a definite 
trend in the variation of pitching moment coefficient is not seen. 

At present, the theoretical prediction of the phenomena of partial 
leading-edge vortex separation is not possible. The extent of the 
separation has been knovm to depend on the leading-edge geometry, the 
Reynolds number and the wing sweep angle (ref. 1). When the degree of 
partial separation can be predicted, the present method can be used to 
calculate the corresponding aerodynamic characteristics. 


h. Conclusions and Recomp^endations 


A theoretical method has been developed for predicting the 
aerodynamic characteristics of low aspect-ratio wings with partial 
leading-edge separation. The present method has been shown to work 
satisfactorily for cases with complete leading-edge vertex separation, 
where the leading-edge suction is zero* Some preliminary theoretical results for 
cases with partial leading-edge vortex separation appear tc be 
reasonable. The method has an advantage over all previous vortex 
lattice methods in that the leading-edge boundary condition con be 
exactly satisfied. It is not restricted to incompressible flew. At 
the present time it is restricted to planforms with pointed wing tips 
only. The present method can be extended to handle arbitrary planforms 
at high angles of attack, as long as the vortex bursting does not 
occur . 

The recommended topics of further research on this methed ore: 

a. Extend the methed to include the wing-tip vertices. 

b. Search for a better iteration scheme for faster convergence 
and lock for a better convergence criteria. 

c. Modify the method for thick wings, 

d. The method can easily be extended to complex planforms in 
which the inboard portion has separated flow and the 
outboard has attached flow. 

e. The method should be checked for some more cases with partial 
leading-edge separation , 
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f . The computer program coding should be made more 
efficient . 
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6. Appendices 

6.1 Appendix A: Evaluation of Induced Velocity Due to a Line Vortex- 


Segment 

In the linearized compressible flow the velocity field induced by 
a line vortex segment of strength T (figure 12) is given by (ref. 31), 


V(R) 


- R) X d£ 
Utt S, ~3 




where 




R = xi + yj + zk 


= Ci + nj + ;k 


R* = xi + gyj + gzk 


= 5i + 6nj + Beit 

Rg = K^-x)^ + B^(n-y)^ + B^(?-z)^' 



R' 


ihe substitution R. - R - a + t£., reduces equation (A.l) to, 


i+2 


U3 


V(R) = axis 

4TT o ,-_2 


dx 


(at"" + Bt + 


6^r 7 r 2 B 

^ (B^ - UA 


2(2 A + B) 


(B^ - 1*AC)(A + B + 


■}, B^ - UAC 0 


(A.2) 


where A = |i,’|^, B = 2a' • 1' and C = |a' |^. Further, it can be 
shown that (ref. 31), 


B^ - 1*AC = -U|a' X t,' 

(A. 3) 

2A + B = 2 b' ♦ 1' 

(A.M 

,-► ,2 

A + B + C = |b' 1 

(A. 5) 


where 

a = (x^ - x)i + (y^ - y)j + - z)k 


b = (x^ - x)i + iy^ - y)J ■*■ (^2 " 


t = (x^ - x^)i + (y 2 - 


a' = {x^ - x)i + 6(y, - y)j + 3(zj^ - z)l^ 
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uu 


b' = (x^ - x)i + 3(72 " 


’ = (x^ - x^)i + 6(72 ■ ^(^2 ~ 


On rearranging, equation (A. 2) becomes 


V(R) 






t' 


(A. 6) 



6.2 Appendix 3: Derivation of Eicpresalons for Press\ire Distribution 

Consider two adjoining sets of spanwise strips of bound elements 


(Sketch B.l). Along the common edge, there are three trailing 
vortices: one due to right set of bound elements, one due to the 

left set of bound elements and the other due to the leading-edge vortex 
system. The force acting on the chordwise element of length Ax of 
the leading— edge vortex system is, by the Kutta— Joukowski theorem, 

F. = pV^r.v.Ax (B.l) 

“ J 1 



i*6 


where p is the fluid density, the free stream velocity, F 

the strength of Jth leading-edge vortex element and v^ the si dewash 
at point i . It follows that the force acting at ith point per unit 
dynamic pressure and length is, 



(B.2) 


where 1 ^ PV^, the diy-namic pressure. A similar expression for the 

force per unit dynamic pressure and length can be written for the 
outside leg of the jth strip as. 



-2v. / 
1 


Y dx 


(B.3) 


where Y is the bound vortex density the leading— edge 

x-coordinate of the trailing-leg under consideration. The 
transformation , 


c . 

X = X. + (1 - Cos 6) 


(B.4J 


reduces equation (B.3) to the form. 


'N. ' 

i_ 

qAx; 


e. 

- ^ ^ Y Sin 6 d6 


IT c . V . 




i-l 

I Yu Sin 8 + 

k=l ^ 




Y- Sin 6. 


:b.5) 



where c is the local chord, N the number of bound elements in 


chordvise direction and 0. 


= - M ^ 


The integral has been 


k 2N 

reduced to a finite s\iri through the regular trapezoidal rule. 
Similarly, for the left leg of (j+l)th strip, 


Wl ' “ 


r 

i-l y. Sin 9. 

\ \ * ' 2 ' 


(B.6) 


Therefore, the force per unit dynamic pressure and per unit length at 
the ith point is given by the sura of equations (B.2), (B.5) and (B.6), 



(B.T) 


Equation (3,7) is evaluated at all endpoints of wing bound 

elements and linear interpolation is performed to obtain the force 

acting at the control station which is inside the vortex strip. Let 

it be denoted by H. . for ith bound element of Jth strip. Then, 

J » ^ 

the contribution to differential pressure coefficient, 
chordvise vortices is, 



T 



(B.8) 


where hy^ is the width of Jth spanwise strip. 

Contribution to AC due to bound elements is calc'olated in the 

P 

following manner. The normal force per unit length acting at ith 





1*8 


bound element of Jth strip (Sketch B.2) is. 


or 


2 

■'b. . ' ♦! - »i> c^TiT 

J .1 1 


Ffi = 2q(u. - V. tan *;^^)Y^Ay 

j .i ^ 


(B.9) 


where is the bound vortex density, and the x and y 

components of the velocity, the sweep angle of the bound elements 

and Ay, the width of the jth strip. It follows that the AC due 
J P 

to the ith bound element is given by dividing equation (B.9) by 

( qAy^ ) ; 



V. tan C. )y . 
1 ' 1 


(3.10) 


The total AC 

P 


is given as the sum equations (B.6) and (B.IO): 



(B.ll) 


Up to this point AC^ has been calculated at the regular wing 

P 

vortex locations. The contribution from the leading-edge vortex 
element on the planform near the leading-edge (EF in figiire 11) has 
yet to be considered. This is done in two steps; 1. Extrapolate AC 

P 

due to the wing vortex system to obtain the AC^ at the location of 
the leading-edge vortex element EF ; and 2, Subtract AC^ induced cy 
the leading-edge vortex element. 



U9 



Sketch B.2 


To obtain AC at any chordvise location, AC sin 6 will be 
p P 

Fourier-analyzed. The factor, Sin 0, is included to eliminate the 
known square root singularity of AC^ at the leading and trailing 
edges. Therefore, let 

N 

AC Sin 6 = a + I a, Cos A0 (B.12) 

P o . I 


where , 


a = — AC Sin 6 d6 
o TT o p 


1 

N 


N 

Z 

k=l 



Sin 6^ 
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a„ = - AC Sin 9 Cos J19 d6 

Jc 7T O P 


^ 2 


N 


rr Z AC Sin 6, Cos £6, 

» Pi, k k 




> k — 1,2,- 


U - Nunter of chordvise lines 

The integrals for Fourier coefficients are reduced to finite sues 
through the mid-point trapezoidal rule. Equation (B.12) can nov be 
used to cs-LCulate the location of the leading— edge vortex 

element EF , which is located at 6 = tt/ 2(N + 1;. To achieve-the 
second step mentioned above, the constant vorticity of the leading- 
edge vortex element is first converted tc vortex density. The 
concentrated vorticity is related to vortex density by, 


r = / Y dx 


(B.13) 


On using equation (B.M in the above equation. 


c . 

r, = -^ / Y Sin e de 

J ^ 


(B.lh) 


Assuming that the concentrated vorticity due to leading-edge vortex 
system is distributed near the leading— edge only and using the mid-point 
trapezoidal rule, the equation (B.li^) reduces to: 



51 


or 


TTC , 


■■j ' 2(ir^ ''j "i 


^ 2(N + 1) j, 

” TT Cj Sin 6^ J 


(B.15) 


where 0^^ = tt/ 2 (N + l). Therefore, the decrease in AC^ value at the 
leading-edge vortex element is given by using equation (B.IO) as. 


(AC ) = - (u - V tan ip)- y. (3.16) 

^ decrease 

where the subscript 2-e means that the variables u, v and are 

evaluated at the leading-edge vortex element. Note that this decrease 

in AC value near the leading-edge from the usual AC distribution 
P P 

is a result of the leading-edge Kutta condition. Hence, the equations 
(B.12) and (B.ll) can be used to calciilate the actual AC^ at the 
location of the leading-edge vortex elements. 
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Number of spanwise strips, (M— l) 
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lU. Variation of number of spanwise strips with aspect ratio 
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Figure 17. Variation of lifo and pitcl.ing moment (about C...5C) coefficients 
with angle of attack for aspect ratio I.II 7 del :.a wing 
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Figure 19b. Variation of pitching moment (about apex) coefficient with 
angle of attack for aspect ratio 2.0 delta wing 










Fig-are 20b. VcLria'iion of pitching moment (about apex) coefficient vrith 

angle of attack for aspect ratio 1.5 delta ving at 0. and 0 , 
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Present method 
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Figtire 22c. AC^ distribution for aspect ratio l.lli? delta ving 
at 30.7 degree suigle of attack 
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Figure 23b. AC distribution for asoect ratio 1.14559 delta vine 

p 

at 19.1 degree angle of attack 
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Figure 23c. 


2y/b ' X ) 

uCp iistriOutio*. for aspect ratio 1 
at 23.9 degree angle of attack 




1.0 


iilia iltr i.: 


Potential (ref. 32) *^v'. 

Present method, 50/° suction lost • 

Present method, 15 % suction lost 
Present method, 90 % suction lost 
Present method, 100% suction lost / 

Suction analogy (ref. 32)^. 

100% suction lost 




a (degrees) 

Fiff.a-e 2i.a. Variation of lift and pitching moment (about 0.25C) coefficient 
with angle of attack for aspect ratio 2.0 delta wing for 
(ixffprent oi xcading-edge suction 
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Present method, 5055 suction lost 
Present method, 75/^ suction lost 
Present method, 90% suction lost 
Present method, 10055 suction lost 
Suction analogy (ref, 32 ) 

100,^! suction lost 
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\^ariaticr. of induced drag coefficient with angle cf atT. 
for aopf.C' ra^io dexta wing for different ariounts o 
leaaing-ea.-e suction 






1. Report No. ^ Government Accession No. 

NASA CR- 145304 

3. Recipient's Catalog No. 

.Viitiial^^^bt.tie A THEORETICAL INVESTIGATION OF THE 
AERODYNAf^lCS OF LOW-ASPECT-RATIO WINGS WITH PARTIAL 
LEADING- EDGE SEPARATION 

5. Report Date 

January 1978 

6 Per^ofrrung Organization Code 

7 Author(s) 

Sudhir Chandra Mehrotra and C. Edward Lan 

8 Perfof ming Organization Report No 

CRlNC-FRL-266-1 

10 Work Unit No. 

505-11-16-07 

fhe*^lyhTver'^iy*’o’i^Taii« for Research, Inc. 

2291 Irving Hill Road - Campus West 
Lawrence, Kansas 66044 

1 1 Contract or Grant No 

NSG-1046 

13 Type of Report and Period Covered 

[ Contractor Report 

12 Sponsoring Agency Name and Address 

National Aeronautics and Space Administration 
Washington, DC 20546 

14 Sponsoring Agency Code 

— — 1 

1 5 Supplerneniary Notes 

Langley technical monitor - Dr. John E. Lamar 


A numerical method is developed to predict distributed and total aerodynamic 
characteristics for low aspect-ratio wings with partial leading-edge separation. 

The flow is assumed to be steady and inviscid. The wing boundary condition is 
formulated by the Quasi-Vortex-Lattice method. The leading-edge separated vortices 
are represented by discrete free vortex elements which are aligned with the local 
velocity vector at mid-points to satisfy the force free condition. The wake behind 
the trai ling-edge is also force free. The flow tangency boundary condition is 
satisfied on the wing, including the leading- and trai ling-edges. Comparison of the 
■*' predicted results with complete leading-edge separation has shown reasonab ly good 
agreement. For cases with partial leading-edge separation, the lift is found to be 
highly nonlinear with angle of attack. 


17 Key Words (Suggested by AuthorUM 

Leading-Edge Separation 
Vortex Elements 
Leading-Edge Thrust 
Low Aspect-Ratio Wings 
High Angle- 0 f-Attack 

18. Distribution Statement 

Unclassified - Limited 
Star Category - 02 

19. Security Oaoit (of this report) 
Unclassified 

70. Security Classif (o< this page) 

Unclassified 

? 1 No of Pages 

83 

77 Pf'ce* 

$6.00 


* For sale by (he National Technical Intoinvation Service. Springfield Virginia 22161 







